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Abstract 

The spiral core instability, observed in large aspect ratio Rayleigh-Benard 

convection, is studied numerically in the framework of the Swift-Hohenberg 

equation coupled to a large-scale flow. It is shown that the instability leads 

to non-trivial core dynamics and is driven by the self-generated vorticity. 

Moreover, the recently reported transition from spirals to hexagons near the 

core is shown to occur only in the presence of a non-variational nonlinearity, 

and is triggered by the spiral core instability. Qualitative agreement between 

the simulations and the experiments is demonstrated. 
PACS: 47.54+r, 47.20.Hw, 05.40.+j, 47.27.Te 
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One of the most intriguing and unexpected recent discoveries in natural pattern formation 
is the observation of spatio-temporally disordered spiral and target patterns (l],|2j in large 
aspect ratio Rayleigh-Benard convection (RBC) in Boussinesq fluids, in a parameter range 
where previously only rolls were known to be stable ||. This regime is characterized by the 
spontaneous and continuous emergence and annihilation of large extended spiral and target 
patterns. Theoretically, these novel states were successfully reproduced both by numerical 
simulations of the Swift-Hohenberg (SH) model coupled to a self-consistent large-scale flow 
(see Eqs. [TJ-^J) |3]|| , as well as by the integration of the full thermally driven Navier-Stokes 
equations in the Boussinesq approximation ||. It is currently postulated that the large- 
scale flow is necessary for the spatio-temporal chaotic state with many spirals and targets 
0-0]. In a first attempt to understand these states, Cross and Tu proposed a physical 
mechanism based on wavenumber frustration, in which defects have an invasive nature and 
create spirals and targets 0. In detailed experiments by Assenheimer and Steinberg |||9|], 
and more recently by Plapp and Bodenschatz |H]] , a new instability of spiral cores of single- 
and multi-armed spirals was observed. The striking feature of this instability is that spiral 
cores oscillate periodically with a frequency considerably higher than the frequency of the 
overall spiral rotation jlUfl . For yet higher supercriticality, a novel transition from spirals to 
hexagons was found in which both up- and downflow hexagons coexist HH . These hexagons 
often invade the background RBC pattern, originating mainly from extended pattern cores. 

In this Letter we present numerical simulations of the spiral core dynamics performed 
in the framework of the SH model coupled to a self-consistent large-scale flow. We show 
that this simple model exhibits both the spiral core instability and the spiral-to-hexagon 
transition, and that both are linked via the large-scale flow. In the case of the spiral core 
instability, we demonstrate that the velocity field generated by the spiral tip increases the 
local wavenumber and eventually drives the tip into the skew- varicose unstable region. Phase 
slips then occur, locally unwinding the spiral and returning the wavenumber into the stable 
domain. This instability is found to have a well-defined threshold in Rayleigh and Prandtl 
number space. Furthermore, the interpretation of the numerical results is supported by 
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the analysis of a similar but simpler problem: a single-armed spiral in an external velocity 
field created by a point vortex located at the spiral core, in the limit of infinite Prandtl 
number when the large-scale flow and the order parameter are decoupled. In addition, 
we observed, at higher supercriticality and only in the presence of non-variational nonlinear 
terms, a transition from spirals to hexagons. Both up- and downflow hexagons are generated 
simultaneously near the spiral core. In this scenario the core oscillations precede and trigger 
the transition to the hexagonal state. 

We considered the well-established model, which describes RBC in a Boussinesq fluid 
rather successfully |TT| , 

$t + (u • V)^ = eV> - 9^ + 3(1 - g){V^) 2 W 2 ^ - 

-(1 + V 2 ) 2 ^ (1) 

fit - a(V 2 - c 2 )Vl = g m z • V(V 2 V>) x Vip (2) 

Q = V x u (3) 

Here ip is the order parameter, u the horizontal velocity field of the large-scale flow, and 
Q the vertical component of the vorticity. The control parameter e represents the reduced 
Rayleigh number, while a characterizes the Prandtl number of the fluid. The parameter g 
allows to more accurately reproduce the stability properties of convection patterns, while g m 
characterizes the coupling strength between the order parameter ip and the vorticity Q. The 
phenomenological parameter c is introduced to describe the local dissipation of the vorticity 
(e.g. due to friction at the bottom of the convection cell) |IT|J12]] . Thus, Eq. || describes the 



dynamics of the order parameter ip, while Eq. [|, using the definition of the vorticity (Eq. 
represents the coupling of the large-scale flow field u and the order parameter. For g = 1 
and g m = Eqs. reduce to the Swift-Hohenberg equation (SHE). 

We solved Eqs. in a domain of 256 x 256 mesh points using a pseudo-spectral method 
based on the Fast Fourier Transform. The physical domain size was typically restricted to 
150 x 150. Circular boundary conditions were enforced by ramping e towards negative values 
at distances r > R ma x = 55. The computations were performed on a parallel Cray J932 



supercomputer, and verified on a 512 x 512 grid. 

We started the simulations from initial conditions of the form x/j = cos(gr + nO), where r 
and 9 are the polar coordinates, q is the wavenumber, and n = ±1, ±2, ... is the topological 
charge of the spiral (|n| is the number of spiral arms, while the sign corresponds to the 
chirality). These initial conditions relax in about 10 to 20 horizontal diffusion times to 
spirals. For sufficiently small values of the parameters e and g m (see Fig. [l|), the spirals 
maintain a stable rigid rotation with an angular velocity depending both on g m and e, as 
well as on the topological charge n ||13|| . Typical spatial distributions of the order parameter 
and the vorticity field for one-armed spirals are shown, e.g. in Ref. ||. As can be seen, 
the spiral tip generates a highly localized vorticity peak at the core. In the case of n-armed 
spirals, n identical vortices are created at the core. In this regime spiral cores are stable and 
only experience a slow off-center drift if the aspect ratio is not sufficiently large Jl2| . 



For e above some threshold, depending on g m (see Fig. p, we observed a novel spiral core 
instability. In contrast to the off-center drift, it persists even in an infinite system, since 
the unstable mode is localized near the core. The main feature of the instability is that 
the spiral core oscillates in the reference frame of the rotating spiral. The amplitude of the 
modulation, caused by the core oscillation, increases from zero as the bifurcaton is crossed 
towards higher values of e and g m . The critical value of the control parameter e depends on 
g m as well as on n. Figure [l] shows the bifurcation diagram as a function of e and g m for 



n = 1,2 [13]. Figure presents typical numerical as well as experimental snapshots of the 



core oscillation for one-armed spirals []14j. Similar behavior occurs for two-armed spirals. 



The core oscillations can be illustrated by comparing the temporal behavior of the order 
parameter at points near the core and at the rim. At the edge, the core oscillations are 
negligible so that mainly the background rotation is sensed. Figure |3| shows plots for one- 
armed spirals for three values of e. Figure |3|(a) illustrates the rigid spiral rotation, below the 
threshold for core oscillations. The central (solid line) and peripheral points (dashed line) 
oscillate with the same frequency, the unambiguous signature of the spiral's rigid rotation. 
Figures |](b) and ||(c) show the spiral dynamics above criticality. Here, the core oscillates 



at a much higher frequency than the peripheral point. The core thus has a fast rotation in 
the framework of the spiral's overall rotation. 

The corresponding experimental data are shown in Fig. [|. Figure f|(a) shows the rigid 
body rotation of a one-armed spiral below the onset of the core instability. Clearly, the 
oscillations near and away from the core are phase-locked. Figure |](b) presents the fast 
oscillations of the core above the threshold of the instability. On this short time scale 
the peripheral signal does not vary significantly and is not shown. Similar dynamics were 



reported in Refs. pJlO 



Our numerical simulations of the full model (Eqs. |TJ-^J) suggest that the vorticity generated 
by the spiral tip plays a major role in the core dynamics. Therefore, the origin of the 
core oscillations can easily be understood - at least qualitatively - in the framework of the 
following simplified model. Consider a one-armed spiral solution of the form A(r) cos (ft, 
where (ft = qr + 9, in the framework of the SHE (Eq. p, coupled to an external velocity 
field generated by a fixed point vortex with circulation T, placed at the center of rotation. 
The sign of the vortex is chosen to correspond to the self-generated vortex in the full model. 
Using the phase approximation (i.e. neglecting the variations of the amplitude A(r)), we 
obtain (see Refs. |]T!?| , |r5f ), that (ft t = T/r 2 do(ft + .... Therefore, the local wavenumber 



increases linearly, q t oc T/r 2 . In other words, the external velocity field winds the spiral 
up near the core. Eventually, the local wavenumber will be carried away from the stable 
band and the skew-varicose (or Eckhaus) instability will be initiated. This stage can no 
longer be described within the phase approximation. Abrupt phase jumps by 2tc (phase 
slips) consequently emerge and return the wavenumber back into the stable region. This 
process then recurs, leading to quasi-periodic oscillations. 

We simulated the spiral dynamics using this simplified model, with a velocity profile 
of the form u = Tr -1 ^, and observed the abovementioned scenario. As the magnitude of 
the circulation T increases, a bifurcation similar to the one described above occurs. At 
small T a steady rotation persists, while for T > r c the core starts to oscillate. Interesting 
dynamics also develop if the sign of the external circulation and the spiral-generated vortex 



are opposite. Then, the local velocity field unwinds the spiral, driving it into the zig-zag 
unstable region fllEf . Because it is non-local, this instability cannot be saturated by phase 
slips and eventually destroys the spiral. This might occur when a spiral is driven by the 
external vorticity created by another oppositely charged spiral and could be relevant for the 
persistence of the spiral turbulent state [|IJ,0]. 

Recently Assenheimer and Steinberg reported a transition from the spiral and target 
chaotic state to a state of up- and downflow hexagons, as the supecriticality e increased 
||. These hexagons started to develop and invade the system primarily from spiral and 
target cores and other defects. Dewel et al. |l7j demonstrated that in the framework of 



the SHE at large e, the coexistence and linear stability of up- and downflow hexagons is 
caused by the excitation of a quasineutral zero mode, which breaks the inversion symmetry 
ip — > —ip. Despite their linear stability, these hexagons are nonlinearly unstable because 
their free energy, in the framework of the pure variational SHE [g = 1, g m = 0), is higher 
than that of rolls. As a result, nuclei of hexagons, immersed in rolls, ultimately shrink. 
Because the generalized SHE (Eq. |IJ-[3|) with either g ^ 1 and/or g m ^ is non- variational, 
a simple relative stability analysis of rolls versus hexagons becomes impossible. 

Numerics performed with a non-variational coefficient 1 — g — 0.25, relatively large cou- 
pling to the vorticity g m = 10 (customary for this model, e.g. Ref. fl2"l ) and supercriticality 
above the threshold for the spiral core instability, e = 1.9, show that hexagons indeed invade 
rolls. We have observed the simultaneous nucleation of up- and downflow hexagons at the 
core which subsequently spread out (see Fig. |5|). Similar dynamics, obtained experimentally, 
is shown in Fig. ^| However, one might speculate that a local wavenumber change near the 
core, rather than the large-scale flow, causes the transition. To compare, we performed 
simulations for the same parameters but without a large-scale flow (i.e. g = 0.75, g m = 
and e = 1.9). In this case spots emerge from the side wall rather than from the core region. 
It is thus plausible to suggest that the large-scale flow, arising primarily near spiral cores 
and other defects, generates the zero mode of the order parameter [JT8J. Thus, the spiral 



core instability initiates the nucleation of hexagons. Since the ip — *> —ip symmetry is still 
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preserved, both up- and downflow hexagons nucleate simultaneously. 

Summarizing, we studied effects of large-scale flow on the dynamics of a spiral core. Al- 
though our computations were performed in the framework of a simplified phenomenological 
model, two characteristic features of the dynamics were found also observed experimentally: 
the spiral core instability and the spiral-to-hexagon transition. The vorticity field generated 
at a spiral core when g m ^ 0, plays a major role in the spiral oscillations. These oscillations 
are not observed in the variational model in which the coupling with the vorticity mode is 
absent (i.e. g m — 0). On the other hand, very similar oscillations are observed in a model 
with a fixed vortex pinned at the center of rotation of the spiral. At higher supercritical- 
ity the vorticity drives the zero-mode of the order parameter ip and prompts the hexagon 
formation. In that case, the core oscillations initiate the transition to the hexagonal state. 
Certainly, a more detailed analysis of a physically more justified (but more complicated) 
model, similar to that of Ref. ||, is desirable. 
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FIGURES 

FIG. 1. Stability diagram for one- (a) and two-armed (b) spirals for c 2 = 2, a = 1 and g = 1. 
As e increases, the spiral core becomes unstable at the dashed line, regaining stability at the solid 
line as e decreases. 

FIG. 2. Snapshots of periodic spiral core oscillations. Top row: simulations with 
e = 1.45,3 = 0-9) 9m = 27, c 2 = 2 and a = 1; time delay between frames 5, and integration 
domain radius R = 55. Bottom row: experiments with e = 2.88 and a = 4.5. 

FIG. 3. Order parameter if; at the core (solid) and periphery (dashed) of a one-armed spiral 
with g = l,g m = 95, a = 5 and c 2 = 1: (a) below threshold (e = 1.1); (b),(c) above threshold 
(e = 1.3 and 1.35). 

FIG. 4. Shadowgraph intensity at the core (solid) and periphery (dashed) of a one-armed 
spiral: (a) below threshold (e = 2.05, a = 4.5) and (b) above threshold (e = 2.88, a = 4.5). In (b) 
only the fast core signal is shown. 

FIG. 5. Hexagon invasion from a spiral core obtained from the full model with 
e = 1.9, g = 0.75, g m = 10, c 2 = 2 and a = 1. Snapshots taken at t = 10, 110, 470, 650, 340, 2350. 

FIG. 6. Experimental hexagon nucleation at a spiral core with e = 3.19, a = 4.5 and time 
delay between frames At = 3.6, 3.6, 22.7, 18.0, 10.7 t v . 
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